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ABSTRACT 

o , 

rvq I For collections of particles in a thermal bath interacting with an asymmetric substrate, it is possible for a 

ratchet effect to occur where the particles undergo a net dc motion in response to an ac forcing. Ratchet 

O I effects have been demonstrated in a variety of systems including colloids as well as magnetic vortices in type-II 

s^ ■ superconductors. Here we examine the case of active matter or self-driven particles interacting with asymmetric 

I — I substrates. Active matter systems include self-motile colloidal particles undergoing catalysis, swimming bacteria, 

artificial swimmers, crawling cells, and motor proteins. We show that a ratchet effect can arise in this type of 

_^ ^ system even in the absence of ac forcing. The directed motion occurs for certain particle-substrate interaction 

■ fj . rules and its magnitude depends on the amount of time the particles spend swimming in one direction before 

1^ ' turning and swimming in a new direction. For strictly Brownian particles there is no ratchet effect. If the 

S . particles reflect off the barriers or scatter from the barriers according to Snell's law there is no ratchet effect; 

however, if the particles can align with the barriers or move along the barriers, directed motion arises. We 

C^ . also find that under certain motion rules, particles accumulate along the walls of the container in agreement 

0^ ' with experiment. We also examine pattern formation for synchronized particle motion. We discuss possible 

applications of this system for self-assembly, extracting work, and sorting as well as future directions such as 

2 ' considering collective interactions and flocking models. 

Keywords: Colloid, optical traps, active matter 

Q ; 1. INTRODUCTION 

There have been a number of experiments examining the behavior of colloids interacting with ordered substrates 
such as two-dimensional square and triangular trap arrays created by optical means.'^''^ Many of the observed 
orderings for repulsively interacting colloids on periodic substrates have similar structures to those of vortices 
in type-II superconductors with periodic arrays of pinning sites.'^ni^J xhe colloids and vortices can be driven 
over the substrates via an applied current, fluid flow, electric fields, or other optical means.'^^'^Il n jg also 
possible to create systems with asymmetric substrates which produce directed motion of particles even in the 
absence of an external dc drive.1^1211 gy applying an ac drive or by flashing the substrate, a net dc flow of 
f~^ ■ particles can be induced which is termed a ratchet eSectP^ In two dimensions ratchet effects can occur in the 

^^ \ presence of symmetric substrates when some other symmetry is broken, such as with a rotating ac drive^SM^Sl 

or by dynamically switching the substrate in a certain order.'^SEU The ratchet effect has been demonstrated 
for colloids on asymmetric substrates^^^ ^"^ and symmetric substrates with asymmetric drives ,'2S1EI1 for cold atom 
systems with flashing substrates ,I32IE31 for vortices in type-II superconductors with asymmetric pinning sites,l2St2ll 
for vibrated granular matter on asymmetric substrates' ^^ " '^°' and for driven liquid drops on asymmetric heated 
plates.!^ Applications of the ratchet effect include sorting of different particles which ratchet at different rates 
C^ ' over the substrate. In the overdamped single particle limit the ratchet effect generally occurs in a single direction; 

however, in collections of interacting particles, ratchet reversals can arise where the dc flow reverses as a parameter 
is varied.lS^IMIIlolllIl 

In all the ratchet systems mentioned above there is some form of external drive, such as a single ac drive, 
flashing of the substrate, or multiple external oscillating drives. There are many systems, termed active matter 
systems, that exhibit motion in the absence of any external drive. Examples include self-catalyzing colloidal 
particles,'*^ microscopic artificial swimmers,'^ swimming bacteria, and active cells.HSl Recently, experiments were 
conducted for swimming bacteria in the presence of an array of asymmetric funnels.'^ No directed motion occurs 
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for Brownian particles or nonswimming bacteria, but when swimming bacteria are added, a ratchet effect occurs 
and the bacteria concentrate on one side of the funnel array. Due to the complexity of the bacterial system, there 
are a wide variety of possible mechanisms for the concentration including hydrodynamic interactions, chemotaxis, 
or collective motions of the bacteria; however, a simple model of individual particles undergoing run and tumble 
dynamics along with a wall- following rule reproduces all of the experimental observations including the buildup 
of particles in the funnel tips and along the wallsPS In this model, the particles move in a single direction over 
a run length Ir during a fixed run time t^, which is the time between tumbling events. Every r^ simulation 
time steps, the particles randomly reorient and move in a new single direction over a distance Ir until the next 
tumbling event a time r^ later. When the particle encounters a barrier, it moves along the barrier with a speed 
determined by the component of its running motion that is parallel to the barrier until the next tumbling event or 
until reaching the end of the barrier. If t^ is small, Ir is small and the particles act more Brownian-like, with no 
rectification occurring for the limit of infinitesimal t^. The addition of steric interactions or thermal fluctuations 
reduces the rectification, while the rectification is enhanced when Ir is increased. With a similar model it was 
numerically shown and demonstrated in experiment that an asymmetric fly-wheel rotates in a preferred direction 
in the presence of a bacterial bath.'*'^^^'^ Other theoretical and numerical worU^has demonstrated rectification in 
swimming bacteria that obey run-and-tumble dynamics under conditions in which, due to long running lengths, 
the particles accumulate and run along the walls. This work assumed no hydrodynamic interactions indicating 
that hydrodynamic effects are not required to produce rectification. The authors also mention that for bacteria 
that collide elastically with the walls, no rectification occurs. Other experiments have also shown directed cellular 
motion over asymmetric substratea^^HSll and even different directions of motion for different cell types on the 
same substrate.'^ Unlike the bacterial case, the mechanics of how the eukaryotic cells move is very important 
and the latter cannot be modeled simply by point particles. These studies open a new field of controlling active 
matter using periodic and asymmetric substrates. If the motion of self driven particles can be well controlled, it 
may be possible to use the particles to perform microscopic work such as transporting larger obstacles or cargo. 
Additionally, biological sorting could be achieved via ratchet effects. In this work we investigate the rules for 
run-and-tumble particles which generate a net dc motion of active matter in the presence of asymmetric barriers. 
We find that if the particles are fully aligned with the barriers, a much stronger rectification effect occurs. If the 
particles refiect or scatter off the barriers than no ratchet effect occurs even for very large Ir ■ This shows that 
the rectification requires the breaking of detailed balance in the interactions with the asymmetric barriers. We 
also find that there is a build up of particles along the barriers similar to what is observed in experiment .33 This 
build up is due to the wall following rule combined with the finite size of the system. If the tumbling times of 
the particles are synchronized, we find that patterned moving density fronts of particles occur after the particles 
accumulate in corners of the containers or in the funnel tips. 

2. COMPUTATIONAL MODEL 

Figure 1 shows an image of our system which consists of an L x L two-dimensional box of length L — 99 with 
confining walls in the x and y directions containing Na active particles. In the center of the system is an array 
of asymmetric funnels, placed in the same geometry used in the experiments for the swimming bacteria.'^ The 
particle density is p = Na/ Lp' ■ The particles move in an overdamped media without hydrodynamic interactions 
according to the overdamped equation of motion 

r7^ = rr(t) + Ff + Ff + rf. (1) 

Here 77 = 1 is a phenomenological damping term. The time dependent motor force F™(t) is represented by run 
and tumble dynamics. Here |F™| = 2.0 while F™ is selected randomly and changed every r^ simulation time 
steps. The resulting run length Ir — Tr5t\F'^\ where 6t = 0.005 is the simulation time step size. Fig. 2(a) shows a 
single active particle with r^ = 1000 and Ir = 10 moving in straight paths for extended distances before randomly 
orienting and moving in other directions, producing a time average velocity of zero. With smaller r^ the length 
scale over which the motion appears Brownian decreases, as shown in Fig. 2(b) for Ir — 1.0. Many active matter 
systems contain thermal fluctuations which we represent as Langevin kicks with the term F^, where {F^) = 
and {F'I' {t)Fj {t')) = 2rikBTSij5{t — t') where /cb is the Boltzmann constant. Fig. 2(c) illustrates the combination 
of run and tumble dynamics with thermal fluctuations. We initially consider F'^ = 0.0 since in many active 
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Figure 1. A snapshot of the system which consists of a confining box with an array of asymmetric barriers (heavy lines) 
and Na active particles (green dots) that obey different rules for motion and interactions with the walls and barriers. 

matter systems at room temperature, such as bacteria, the thermal fluctuations are small on the scale of the 
particles. The term F^'' is the particle-particle interaction force which is modeled as simple short ranged steric 
repulsion. The term F^ is the force from the barriers which are each modeled by two half-parabolic domes of 
strength Fb — 30 and radius of tb — 0.05 which repel particles along the direction perpendicular to the trap 
axis. 

■i^B \^ ^Bri , ,'± FBr2 , ,' 

^ ^ -e(ri)R,^^ + — ^e(r2)R 



fc=0 



rB 



TB 



'ik ■ 
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Here the total number of barriers including the confining walls is Nb = 28, ri = tb — R^i., ^2 = tb — i?^, 

^ik = \^i -^k ±-^sP||l: R-jfc = (R-i -R-fc ±^-BP||)/^ifc, Rik = l(R-i -R-fc) -Pj^l: R-ifc = [(R-i -Rfc) -vA/Rik- 

Ri(Rf ) is the position of particle i (barrier fc), and pn (p'[) is a unit vector parallel (perpendicular) to the axis 
of barrier k. The funnels are modeled as two barriers meeting at a common end point at angles 9 and t: — 9 with 
the X-axis and the tips of adjacent funnels are spaced a distance 8.25 apart along the x direction. 

3. RULES PRODUCING RECTIFICATION FOR NON-INTERACTING PARTICLES 

We first consider the case of no thermal fluctuations and no steric interactions between particles. In previous 
work, a particle encountering a barrier had the component of its motion perpendicular to the barrier canceled 
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Figure 2. The trajectories of a single active particle undergoing run and tumble dynamics, (a) Ir — 10. (b) Ir = 1. (c) 
The same as in (a) but with an added thermal noise component with F^ = 10.0. 
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Figure 3. Different rules for the particle-barrier interactions, (a) Rule I. The particles move with the component of F"^ 
that is parallel to the barrier, (b) Rule II. F™ is realigned to be completely parallel to the barrier, (c) Rule III. The 
particles reflect off the barriers, (d) Rule IV. The particles scatter off the barriers. 



but experienced no reorientation of F™ by the barrier.HS This behavior, termed Rule I, is is illustrated in the 
schematic in Fig. 3(a). Here we consider three additional barrier interactions. In Rule II, F™ is realigned to 
be parallel with the barrier. This is a realistic assumption for many types of active matter such as elongated 
bacteria, which can align their swimming with a wall due to hydrodynamic interactions. In Rule III, F™ is 
reversed so that the particles reflect back into the direction from which they came, as shown in Fig. 3(c), while 
in Rule IV, F™ is reflected across a line perpendicular to the barrier so that the particles scatter off the barrier, 
as illustrated in Fig. 3(d). 

We now consider Rule I, partial alignment with the barriers, for short and long l^ as illustrated in Fig. 4. 
For Ir = 0.01 after 100 tumbling events in Fig. 4(a), there is no rectification and the density p^^^ j p in the upper 
channel equals the density p^^^ j p in the lower channel, as shown in Fig. 5. For Ir = 180, there is a buildup of 
particles in the upper chamber as shown in Fig. 4(b), and the dependence of p^^Vp ^'^^ P^"^^ I P on time has a 
shape very similar to the form found in experiments, as illustrated in Fig. 5. The rectification mechanism for 
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Figure 4. Image of the sample after 100 tumbling events for Rule I. (a) At Ir = 0.01, the particle density is equal on both 
sides of the chamber, (b) At l^ = 180 the particle density has substantially increased in the upper chamber with particles 
accumulating along the walls and in the funnel tips. 

long Ir is illustrated in Fig. 6. When a particle approaches a funnel from below, it is guided along the barrier and 
moves into the upper chamber as long as r^ is long enough that it does not change direction before reaching the 
upper chamber. On the other hand, when a particle approaches a funnel from above, it is guided into the closed 
tip of the funnel where it remains trapped until a new tumbling event sends it back into the upper chamber. For 
small Ij. the same barrier interaction rules still apply; however, due to the short time between tumbling events, 
the particle spends so little time moving along the barrier that its motion cannot be guided and no rectification 
occurs. This result shows that it is the combination of the long running time and the breaking of detailed balance 
in the particle-barrier interactions that lead to the rectification. Fig. 7 shows a plot of r = p^^^/p*^^-' versus Ir for 
the Rule I sample. The amount of rectification increases with increasing Ir- In the limit of small Ir where the 
interactions between the particles and the barriers are in the Brownian limit there is no rectification. 

When we employ Rule III or Rule IV we find no rectification, as shown in Fig. 8 where we plot r for the 
four different rules in a sample with 4 = 300. Regardless of the value of Ir, Rules III and IV do not produce 
rectification of the particles due to the fact that they preserve detailed balance for the interactions of the particles 
with the barriers, which is enough to prevent rectification. 

We note that the amount of rectification varies with other system parameters such as the angle of the funnels, 
the overall size of the funnels relative to lb, and the spacing between adjacent funnels. In general, if the funnels 
are made larger for fixed Ir the rectification is reduced, since this has the same effect as holding the funnel size 
constant and reducing Ir- Increasing the spacing between the funnels reduces the effect since there is simply 
more space for particles to freely pass from the upper chamber to the lower chamber, producing a larger amount 
of backflow and lowering the efficiency of the ratchet. As a function of the angle between the funnel arms, the 
ratchet effect is generally optimized near 50 to 60 degrees.H3 

4. STERIC AND THERMAL INTERACTIONS 

We next consider the effect of adding a short range steric repulsion between the particles which is modeled as a 
stiff spring of finite radius which produces a force 
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Here r^ — Va — Rij, Rij — |Ri ~ R-jl, R-ij = (R-i — ^j)/Rij, and fa — 150 is the force coefficient. 
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Figure 5. p^^^ I p and p'^'^^ / p vs time in simulation time steps for the system in Fig. 4 under Rule I. Dotted lines: Z,. = 0.01 
from Fig. 4(a) where the density is almost the same in each chamber. Upper heavy line: p I p for the l,. = 180 system 
from Fig. 4(b); lower heavy line: p I p for the system from Fig. 4(b). 





Figure 6. A schematic of the rectification mechanism under Rule I and Rule II. Left panel: A particle approaching the 
funnel from above aligns with the barrier and moves into the funnel tip where it is trapped. After the next tumbling event 
the particle may escape from the funnel tip and move back into the upper chamber. A particle approaching the funnel 
from below aligns with the barrier and is guided into the upper chamber. Right panel: For very short Ir, the particle 
spends little to no time interacting with the barrier and its motion cannot be guided by the spatial asymmetry of the 
barrier. 
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Figure 7. The ratio of the particle densities r — p'^YP "^s Ir for the system with Rule I. For small Zr the rectification 
is lost. 
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Figure 8. The value of r after 6000 tumbling events for a system with Ir — 300 for Rule I (filled circles), Rule II (filled 
squares), Rule III (filled diamonds), and Rule IV (open triangles). Here only Rules I and II produce rectification. 
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Figure 9. (a) r vs p for a system with steric interactions under Rule I at Z^ = 20. Here r decreases with increasing p. (b) 
r vs the thermal force F^ for a system with no steric interactions under Rule I at Ir = 20. The rectification is reduced 
as the magnitude of the thermal fluctuations increases. 



In Fig. 9(a) we plot r vs p for a system with Ir = 20 and steric interactions under Rule 1. Here the rectification 
decreases with increasing density. The interacting particles interfere with each other's motion in several ways. 
Particles that strike the barriers at an angle that is nearly perpendicular to the barrier move very slowly along 
the barrier. These particles get in the way of other particles that strike the barrier at lower angles and that would 
otherwise move relatively quickly along the barrier. This interaction decreases the efficiency of the ratchet. The 
buildup of particles along the walls also becomes much more apparent for the sterically interacting particles, as 
illustrated in Fig. 10. A similar formation of excess density along the walls occurred in the experiments of Ref.!^ 
In Fig. 9(b) we illustrate the effects of adding thermal noise to a system of particles without steric interactions 
under Rule I with Z,. = 20. The plot of r versus F^ shows that the rectification is strongly reduced as the thermal 
fluctuations begin to dominate the motion of the particles and the system enters the Brownian limit. 



5. SYNCHRONIZED SWITCHING AND PATTERN FORMATION 

Up to this point we have considered the effects of random tumbling events. All of the particles had equal 
Tr] however, at the beginning of the simulation each particle was assigned a randomly chosen amount of time 
< t < Tr after which the first tumbling event for that particle would occur. This desynchronized the switching 
of the particles. It could, however, be possible for certain systems to have synchronized tumbling events for all 
particles. For example, this would occur if a periodically applied external field activated the tumbling, so that 
all of the particles would change direction simultaneously. For open systems or samples with periodic boundary 
conditions, pattern formation would not be expected to occur; however, in the confined box and funnel geometry 
that we study we find that pattern formation in the form of expanding fronts of particles occurs when the 
switching events are synchronized for particles that have no steric interactions. In Fig. 10(a,b) we illustrate the 
formation of moving fronts of particles in a system with no steric interactions at Ir = 80 under Rule 1 with 
synchronized tumbling times. Density fronts are emitted by the four corners of the container and to a lesser 
extent also from the ends of the funnels. These fronts move out in a growing circular pattern, decreasing in 
density as they spread until the structures break up due to interactions with the walls. The pattern forms since 
the particles can effectively collapse into the corners and the funnel tips due to their lack of steric interactions; 
the particles are trapped in these locations due to the barrier interaction Rule I. If r^ is long enough essentially 
all the particles will become trapped in the corners and funnel tips. Since the tumbling events are synchronized, 
the particles move out from each trapping point in a circular pattern after a switching event. In our system 
we observe only a semicircle of motion due to the confinement by the walls. The fronts are clearly defined only 
for sufficiently large r^ since for short t^ the particles do not have enough time to become highly concentrated 
in the corners and funnel tips. This type of pattern formation is strongly reduced for particles with steric 




Figure 10. The build up ol particles along the walls and inside the funnels for a system with Rule I and sterically 
interacting particles at Ir — 120. This same effect is observed in experiments. 
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Figure 11. The formation of moving particle fronts for a system with Rule I, no steric interactions, Ir — 80, and 
synchronized tumbling times, shown as consecutive snapshots of the system taken at equal time intervals. The particles 
accumulate in the corners of the container and are released as a semicircular density wave after each switching event. 

interactions since the interactions reduce the ability of the particles to become concentrated at a point; however, 
some remnant of the moving fronts still can be observed when steric interactions are present. 



6. DISCUSSION AND FUTURE DIRECTIONS 

This work in conjunction with experiments shows that simple substrates can be used to create directed motion in 
self-driven particle systems. This opens the possibility for an entirely new class of ratchets based on microscale 
systems which could be used for sorting different species of active particles, guiding particles to move in certain 
directions, and possibly even extracting work out of such systems. Our results show that thermal fluctuations 
can strongly reduce the efficiency of the directed motion. For large colloids and many biological systems, thermal 
effects are generally small on the scale of the particles; however, if nanoscale self-driven particles are considered 
then thermal effects will be relevant. Our work and other models of this type of system include only point particles 
or elongated rods. It would be interesting to study ratcheting behavior for particles with more complicated shapes 
or for particles that contain internal degrees of freedom. For example, eukaryotic cells exhibit different types of 
mechanical mechanisms for locomotion such as pushing or pulling motions, and it could be possible to implement 
swimming rules that mimic such motions. We have studied only static substrates in our system; however, it should 
be possible to create dynamic substrates for self-driven particles similar to flashing ratchets, and such dynamic 
substrates might be much more effective in directing the motion of the particles. It would also be interesting to 
study different types of collective effects such as by introducing models for flocking or swarming behaviors or by 
adding longer range interactions between particles which could arise through hydrodynamic effects. Although 
the model we consider is mostly relevant for microscopic systems, similar effects could be tested for larger scale 
systems such as individual or collectively moving insects, animals, and birds. Another possibility would be to 
study self-driven particles that have an intrinsic asymmetry as they move through symmetric or asymmetric 
substrates. 



7. SUMMARY 

We have shown that a new type of ratchet system can arise in systems of self-driven particles or active particles 
in asymmetric geometries. Unlike most ratchet systems where some form of external ac force or flashing potential 
is required to induce a ratchet effect, the self-driven particle system can exhibit a ratchet effect in the absence 
of external forcing. In our system we consider active particles that move with a run-and-tumble dynamics of the 
type found in bacterial systems where the particles move in a single direction for a fixed time before undergoing 
a tumbling event, randomly reorienting, and then moving for a fixed period of time in a new direction. We 
studied different rules for the interactions between the particles and the barriers, including alignment or partial 
alignment with the barrier which causes the particles to follow the barriers as well as reflection and scattering 
interactions. When the particles are allowed to move along the barriers and when the time between tumbling 
events is sufficiently long, we observe rectification of the particles by the barriers. We also find a buildup of 
particles along the walls and inside the funnel tips. Both the rectification and the buildup of the particles along 
the barriers are in excellent agreement with experimental observations. Under these same wall interaction rules, 
if the time between tumbling events is very small the motion of the particles becomes more Brownian-like on 
the length scale of the barriers and the rectification is lost. For reflection and scattering of the particles off the 
barriers the rectification is lost even for very long times between tumbling events. These results show that it is 
the combination of the running length and the breaking of detailed balance in the particle-barrier interactions 
that produce the rectification. We have also found that if the tumbling times for the particles are synchronized, 
pattern formation occurs in the form of moving density fronts generated after each switching event due to the 
concentration of particles in the corners and funnel tips. The addition of thermal effects or steric interactions 
generally reduces the effectiveness of the rectification. Future directions to pursue include more complicated 
shapes and interactions for the individual particles as well as the addition of collective effects. 

This work was carried out under the auspices of the NNSA of the U.S. Department of Energy at Los Alamos 
National Laboratory under Contract No. DE-AC52-06NA25396. 
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